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Abstract 



We calculate the orbital magnetic response to Aharonov Bohm flux of 
disordered metallic rings with attractive pairing interaction. We consider the 
reduced BCS model, and obtain the result as an expansion of its exact solution 
to first order in the interaction. We emphasize the connection between the 
large magnetic response and the finite occupation of high energy levels in the 
many-body ground state of the ring. 

I. INTRODUCTION 

One of the remarkable phenomena of mesoscopic physics is the existence of equilibrium 
persistent currents in small normal metal rings, in the regime where the elastic mean free 
path / is much smaller than the ring's circumference L. This was both predicted theo- 
retically [0] and observed experimentally. For the ensemble-averaged persistent 
current, the experimental results are much larger than the value obtained using the 
model of noninteracting electrons, and show predominantly diamagnetic response at 
zero flux. In an attempt to resolve this discrepancy, the contribution of electron-electron 
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■'e-e") interactions was calculated |]T0[ and it was later suggested that the diamagnetic re- 



sponse is due to effective attractive e-e interactions. |Tl| Though indeed the inclusion of e-e 
interactions increased the theoretical value, it was still smaller by a factor of about 5 than 



the experimental results, for both repulsive and attractive interactions. [11 



In a recent work [|T^] we suggested that the contribution of high energy levels (denoted 
"far levels"), further than the Thouless energy E'xh from the Fermi energy Ep, results in an 
enhanced orbital magnetic response (the derivative of the persistent current at zero flux). 
This was shown for the model of attractive pairing interactions described by the reduced 
BCS Hamiltonian, by doing perturbation theory in both the magnetic field and the e-e 
interaction. In this paper we repeat the calculation for the same model, doing perturbation 
theory in the magnetic field only, using the exact many-body states of the system with e-e 
interactions. In this way we can treat on the same footing different regimes of the strength 
of the interaction. We consider briefly the limits of zero interaction and the opposite limit 
where the interaction is strong enough and the system is superconducting. We then treat in 
more detail the case of weak e-e interaction, using Richardson's exact solution for the reduced 
BCS Hamiltonian. [0,0 In this way we obtain [see Eq. (|T^)] the result of Ref. [|I^ as a 
leading order expansion of an exact solution, instead of first order perturbation theory in the 
interaction. Furthermore, the present derivation has the merit of emphasizing the connection 
between the enhanced magnetic response and the pairing correlations of all the levels up 
to the high energy cutoff at the Debye frequency cjd. These correlations exist in the exact 
many-body ground state (g.s.) of the ring [see Eqs. ([T^) , (p!5|) and discussion after Eq. (0)]. 
The contribution of the pairing correlations of the far levels was first realized in connection 



to superconductivity in small grains. [|T6[ There it was shown that this contribution results 
in a much larger condensation energy than that given by the BCS theory in a large regime 
in which superconducting correlations are well developed, as well as in a correction to the 
spin magnetization and susceptibility as function of magnetic field H that persists up to 
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II. MAGNETIC RESPONSE 



We consider a quasi one dimensional disordered ring with a finite number of electrons 
A^e; penetrated by a constant Aharonov Bohm (AB) flux in its middle. We take the case of 
zero temperature, and the average spacing between noninteracting energy levels is d. The 
Hamiltonian is given by 



H = Y,j drii(r) 



(1) 



where U{r) is the external potential which includes the disorder, and H^e represents the e-e 
interactions. The vector potential corresponding to the AB flux $ in the middle of the ring 
is given in the London gauge hy A = $/(27rp)0, where p is the distance from the origin and 
the angle (f) is in the clockwise direction of the ring. The ground state energy of the system 
is flux dependent, and can be written for small flux as 



E{^) = Eo- -^2$^ + ... 



(2) 



The persistent current is given by / = —dE/d^. Since, due to time reversal symmetry 
there is no linear term of the energy as function of flux, hence /(O) = and dl/d^ |$=o= E2. 

Since we are interested in E2, we calculate the ground state energy of the system to 
second order in the flux. We take as the unperturbed Hamiltonian 



p2 



V'a(r) + He, 



(3) 



and the magnetic field represented by the vector potential as perturbation 



iP-A + A-P) + 



:A' 



(4) 



2mc^ ' 2mc^ 

We assume that the width of the ring is much smaller than its circumference L. |18 
Then the A"^ term gives the "diamagnetic" contribution. 



1 / 6$ \- 



N 
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which is independent of the e-e interactions and results in 



(5) 



We denote by \i) the eigenstates of the noninteracting electrons in the disordered ring 
without magnetic field. In this basis the first and relevant term of Hi is 

Here q destroys an electron in the state \i) with wavefunction Xi{^) ^"^^ Pij = (^l-^l |j) is the 
matrix element of the momentum parallel to the ring's direction. We choose the Xi's to be 
real, and then Pij is pure imaginary and Pa = 0. Using second order perturbation theory in 
Hj we write the paramagnetic part of E2 as 

ppar e \2 (g.S.|P,,cLc^.Jj) (/ 1 Pfe/cL C,,, |g.S. ) 

I ijkl,aa' 6-°- ^ 

where Ej are the energies of the intermediate states |/). We model the e-e interactions by 
the reduced BCS interaction 



He^ = -\dJ2' 4Al^jiCj^ > (9) 



where A is the dimensionless pairing parameter and the sum is over all levels with energies 
between Ep — ud and Ep + u-p,. This interaction Hamiltonian is the usual one used when 
discussing superconducting grains, both in the perturbative and nonperturbative regimes 
p!9| and its validity is discussed in, e.g. Refs. (In particular, for the model to 

be valid the grain's dimensionless conductance g must be much larger than one.) In this 
model the ground state has no singly occupied noninteracting states |l4| , pr5[ (we assume, for 
simplicity, that the number of electrons in the ring is even). Therefore, |J) has two singly 
occupied states, with opposite spins. We denote by Imn the set of many-body states where 
state m(n) in occupied with one electron with spin up (down). Eq. (P) can then be written 
as 



Er" = -2 



f e Y 
\mcL 



E E \p^n\' ' ; :v ^ • (10) 



We now analyze this equation for the cases of A = (normal metal), A > 1/lnA^ 
(superconductor, see Ref. 0), and < A ^ 1/lnA^ (weak attractive interactions). Here 

For A = the ground state is the noninteracting Fermi state, and the only relevant 
intermediate states are those with one electron-hole pair (the lowest energy states within 
each subspace Imn)- A straightforward calculation results in 

^"■"" = (;li)' S ^ 

where m > denotes states with energies larger than E-p and Umn = Cm — ^n, the differ- 
ence between the energies of the single particle states m and n. For a diffusive ring this 
paramagnetic term is of the same order as the diamagnetic term in Eq. (^. The difference 
between these terms is of the order of the contribution of the last level (and can therefore 
have either sign) and constitutes the noninteracting sample specific result for the magnetic 
response. [|17| 

The opposite limit is the superconducting regime. In this regime one can use the BCS 
approximation for the ground and excited states of the system, and the Bogoliubov trans- 
formation for the creation and annihilation operators. Eq. ( [T(]| ) is then reduced to 

^par(BCS) _ 2 ^ V 'V \^'^^\^ {UmVn — UnVmY /^2^ 

Here Um, Vm are the coherence factors, and Em is the energy of the electron- hole quasiparticle 
of state m. In the ballistic case, Pmn = Pm^mn and therefore the paramagnetic term is zero. 
This results in the well known perfect diamagnetism of a superconductor. For the diffusive 
case the ensemble- averaged momentum matrix elements are given by |13| 



,2\ _ Vl^d 

which is roughly constant for Umn < 1/t and zero for Umn > 1 /t (r is the elastic mean free 
time and s = 1,2,3 is the effective dimension of the ring for diffusive motion). The total 
response in this case (sum of diamagnetic and paramagnetic terms) is diamagnetic, with an 



approximate magnitude of {l/^)Ef^, where ^ is the (balhstic) superconducting coherence 
length. 

We now turn to the calculation of E^"^^ for the case of weak attractive interaction (A < 
1/ In A^). We calculate the interaction correction to E2 to first order in A. Using Richardson's 
exact solution ||TH] one finds that to first order in A the ground state can be written as 



^g.s.= E <P{h...fN)h],...h]^\vac) (14) 

{/i-/iv} 

with 

<j){l...N) = 1 

Here (/i.../Ar) denotes a set of out of the 2N noninteracting eigenstates between Ep —ub 
and Ep + cud, = cj-^tcj-^ and (f){l...N; 7^ j, k) is the amplitude of the many-body state 
with a filled Fermi sea except a pair excitation from state j below Ep to state k above 
Ep. To first order in A, the amplitude of all the other possible configurations is zero for 
the ground state. The finite amphtude to occupy noninteracting states with energies larger 
than Ep is a result of the interaction. For any small A the system gains energy by having 
a different ground state than the Fermi state. Note that the amplitudes of all the different 
configurations come with the same sign. Similar analysis of the intermediate states shows 
that the only intermediate states that contribute to E2 in first order in A are the lowest 
energy states of each subspace Imn, with m and n on opposite sides and within cud of the 
Fermi surface. Then 

(-^mnkmT'^nT ~ '^nl'^mlld ■S-'^ = 1 — — (16) 

and Eg s.—Ej = Umn + ^d. Thus, the pairing interaction contributes to the magnetic response 
in two ways. First, due to the finite occupancy of levels above Ep the matrix element in 
Eq. ( P^ ) has a contribution not only from the annihilation of an electron in state n below 
Ep and the creation of an electron in state m above Ep. Since in the ground state there is a 
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finite amplitude for state m to be doubly occupied and state n to be empty, there is a finite 
contribution to the matrix element from annihilating and electron in state m above Ep and 
creating it in state n below E-p [second term in Eq. (|1^)]. Second, it adds a term Xd to the 
energy of the excited states due to the excess pairing energy of a doubly occupied state. 
The latter term is a Hartree-like contribution, coming from the diagonal part of the pairing 
interaction, while the former is due to pairing correlations of pairs in different single-particle 
states [offdiagonal part of the Hamiltonian in Eq. @]. Both of these contributions suppress 
the paramagnetic term in the case of attractive interaction, and as a result, to first order in 
A, 

Er^'°^ = Er^''^ -xd(^y y (17) 

Using Eq. ([T^) we find that the derivative of the persistent current at zero flux, to first order 
in the interaction, is given by 



^?) = ^lni^. (18) 



In comparison with the known result for the first order interaction correction to the mag- 
netic response [|l^], our result has a much larger logarithmic cutoff, Ub compared with E^h- 
This enhancement is irrespective of the higher order correction, which for attractive inter- 



action further increases the first order result. ||Tl|] The large logarithm we obtain is a result 
of enhanced pairing correlations of all the states within 1/r from Ep. We expect that semi- 
classically this term necessitates only one circulation of the ring, and will therefore affect 
the magnitude of the persistent current as well (see discussion in Ref. [|I|). This may lead 
to a resolution of the discrepancy between the experimental results and theory of persistent 
currents in diffusive normal rings. 
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